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Abstrac^t 

V 

In  this  paper,  we  presentS|,a  unified  method  for  the  delay  analysis 

^  n- 

of  a  large  class  of  random  multiple-access  algorithms.  -0^  method  Is 
based  on  a  powerful  theorem  referring  to  regenerative  processes.  In 
conjuctlon  wlth^esults  from  the  theory  of  Infinite  dimensionality  linear 
systems.  JUe-  apply  the  method  to  analyze  and  compute  the  per  packet 
expected  delays  induced  by  two  algorithms,  in  the  presence  of  the  Poisson 
user  model.  The  considered  algorithms  are:  The  controlled  ALOHA 
algorithm,  and  the  '^0.487'*^  algorithm.  The  same  method  has  been  previously 
applied,  for  the  delay  analysis  of  certain  limited  sensing  random  access 
algorithms .  /c^ 
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1.  INTRODUCTION 

h  key  problem  In  the  design  of  communication  networks  is  the  efficient 
sharing  of  a  common  transmission  channel,  (such  as  a  satellite  link,  a  ground 
radio  channel,  a  computer  bus,  a  coaxial  cable,  or  an  optical  fibre)  among  a 
large  population  of  network  users.  This  problem  is  referred  to  as  the  multiple- 
access  problem,  since  many  independent  users  share,  and,  thus,  access  a  common 
channel  for  transmission  of  information.  The  solution  to  the  multiple-access 
problem  must  Incorporate  a  distributed  control  scheme,  termed  multiple-access 
algorithm,  for  allocating  the  channel  resources  among  the  network  users. 

The  design  and  perfonunce  of  multiple-access  algorithms  are  highly 
dependent  on  the  nature  of  the  users.  VIhen  a  channel  is  to  support  large 
numbers  of  bursty  (low  duty-cycle)  users,  random > multiple-access  algorithms 
(HMkAs)  become  more  efficient  than  d.etern44[i4j3tic_  algorithms,  nils  has  been 
early  recognized  by  the  researchers  in  the  field,  and  a  plethora  of  RMAAs 
have  been  proposed  during  the  past  fifteen  years  [1,2]. 

The  key  performance  measures  of  a  RMAA  are  its  throughput  and  delay 
characteristics.  The  evaluation  of  such  characteristics  has  been  the  subject 
of  numerous  studies.  In  most  cases,  a  Markovian  model  is  employed,  and  the 
existence  of  steady  state  of  the  random-access  system  is  related  to  the 
ergodicity  of  an  underlying  Markov  process.  Depending  on  the  complexity  of 
the  state  space  of  such  a  process,  this  formulation  usually  gives  sufficient 
information  on  the  maximum  input  traffic  rate  that  an  algorithm  can  maintain. 
However,  the  evaluation  of  the  delay  characteristics  is  a  much  harder  problem, 
since  they  are  intimately  itnerwoven  with  the  dynamical  behavior  of  the 
algorithm's  scheduling  mechanisms.  Due  to  this  fact,  it  is  not  surprising  that 
results  concerning  the  delay  characteristics  are  limited,  and  are  obtained  after 
a  rather  Intricate  and  difficult  analysis,  which  is  usually  matched  to  the 
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peculiarities  of  the  specific  algorithm  at  hand. 

In  this  paper,  we  show  how  the  delay  analysis  of  RHAAs  can  be  unified  and 
simplified,  by  the  use  of  ^some  known  results  from  the  theory  of  regenerative 
processes,  and  the  theory  of  infinite  dimensional  systems  of  linear  equations. 
After  outlining  the  method  in  section  2,  we  demonstrate  its  wide  applicability 
and  relative  simplicity,  by  applying  it,  in  sections  3  and  4,  to  two  algorithms 
that  represent  different  classes  of  RMAAs,  namely: 

1)  the  Oontrolled  ALOHA  algorithm  ( "ALOHA-type*  class)  [61 

2)  the  "0.487"  algorithm  ("full  sensing-blocked  access"  class)  [7,8] 

Ebr  the  above  algorithms,  we  obtain  explicit  results  on  the  Induced  mean 
delay,  for  the  Poisson  infinite-user  population  model.  The  higher  moiaents  of 
the  delay,  for  the  Poisson  as  well  as  for  an  arbitrary  memory less  Input  stream, 
can  be  computed  using  the  same  method.  We  note  that  the  method  has  been  already 
applied  for  the  delay  analysis  of  a  class  of  limited  sensing  random  access 
algorithms  111,18,19]. 

2.  THE  METHOD 

In  random-access  systems,  as  in  virtually  every  queueing  system,  many  of 
the  involved  stochastic  processes  are  regenerative .  The  following  definition  is 
taken  from  [4] . 

Definition  The  process  is  said  to  be  regenerative  with  respect  to  the 

renewal  sequence  every  positive  integer  M  and  every  sequence 

^l’*‘’^M’  joint  distribution  of  »..»X^  is  inde- 

11  Mi 

pendent  of  i.  The  random  variables  and  ^i+1  “  referred 

to  as  regeneration  points  and  regeneration  cycles,  respectively. 

For  regenerative  processes,  the  following  elegant  and  powerful  result  holds, 
which  will  be  referred  to  as  the  regeneration  theorem  [3,4,5]. 
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Theorem  1 

Let  the  discrete-time  process  be  regenerative  with  respect 

to  the  renewal  process  Also,  let  -  ^i+i  “  ^i»  1*1 » 2,..., 

denote  the  length  of  the  i-th  regeneration  cycle,  and  let  £  be  a  nonnegative, 
real  valued,  measur2d}le  function. 

If  C  =  bCC^}  <  •  and  S  -  e1  Z  f(Xj|^)}  <  «,  then, 

i=l 

in  1  **  *  S 

fim  -  Z  fCX^)  -  tim  -  B  {  Z  fCX^)  |  »  -  ,  w.p.  1 

n-H>»  n  i«l  n+«»  "  i-1  ^ 

Furthermore,  if,  in  addition  to  the  flnlteness  of  C  and  S,  the  distribution  of 

Cf  is  not  periodic,  then  Xj,  converges  in  distribution  to  a  random  variable 
X,a,  and 

E(f(X«)>  -  - 
C 

Thus,  under  the  conditions  stated  above,  the  limiting  (elected)  average, 
and  the  mean  of  the  limiting  distribution  of  {£(X„)}„>i  exist,  coincide, 
and  are  finite.  Horeover,  their  common  value  is  then  given  in  terms  of  the  per 
cycle  quantities  S  and  C. 

Given  a  RMAA,  let  be  the  process  of  interest  associated  with 

the  random-access  system;  this  process  might,  for  example,  be  the  delay  process 
induced  by  the  algorithm.  Then,  provided  that  {Xn)n>1  c*”  shown  to  be 
regenerative,  the  regeneration  theorem  itself  shows  the  way  to  establish  the 
existence  of  steady  state,  and  to  compute  the  steady-state  moments,  and  the  distribu¬ 
tion  of  (X„}n  by  appropriately  selecting  the  function  f. 

In  virtually  all  existing  RMAAs,  it  is  relatively  easy  to  identify  regenera¬ 
tive  tines  (e.g.,  when  the  system  becomes  empty,  or  when  an  appropriate  Markov 
chain  hits  a  suitable  fixed  state),  at  which  the  process  of  interest  probabilis- 


tlcally  restarts  Itself.  Given  a  RMAA  and  a  function  £,  the  problem  then  is  to 
exploit  the  dynamics  of  the  algorithm,  to  find  those  per  cycle  properties  of  the 
sample  function  of  the  process,  that  could  be  subsequently  used  to  evaluate  the 


quantities  C  and  S. 

In  section  3,  it  is  shown  that  for  the  delay  process,  and  for  f(x)  =  x, 
the  computation  of  S  and  C  are  intimately  related  to  the  solution  of  an  infinite 
dimensional  system  of  linear  equations*  It  can  be  shown  that  this  is  the  case 
when  f(x)  =  x",  n  -  2,3,...,  as  well  [11].  Therefore,  the  steady-state  moments 
of  the  delay  process  induced  by  a  particular  algorithm,  can  be  computed 
from  the  solution  of  the  corresponding  infinite  linear  system.  In  Appendix  A, 
we  give  a  number  of  general  results,  that  are  useful  in  establishing  the 
existence  and  uniqueness  of  a  solution,  and  in  developing  approximations  to 
the  solution  of  such  systems*  In  section  4,  we  apply  these  results  to  the 
specific  infinite  linear  systeois  developed  for  the  two  algorithms  of  section 
3.  This  procedure  involves  the  following  steps. 

Step  1  Find  conditions  under  which  the  infinite  linear  system  has  a 
unique,  nonnegative  solution. 

Step  2  Show  that  the  variables  of  interest  coincide  with  the  unique 
solution. 

Step  3  Develop  arbitrarily  tight  upper  and  lower  bounds  on  the  solution. 

3.  TWO  ALGORITHMS  AND  THEIR  RELATED  SYSTEMS  OF  EgUATIOMS 

For  both  algorithms  of  this  section,  we  assume  a  collision-type,  packet' 
switched,  slotted,  broadcast  channel*  The  channel  is  accessed  by  a  very  large 
(effectively  infinite)  number  of  identical,  independent,  packet-transmitting, 
bursty  users.  The  cumulative  packet  generation  process  is  modelled  as  a  Poisson 
process,  with  intensity  1  packets  per  slot.  However,  the  proposed  method  can  be 
applied  equally  well,  when  the  number  of  packets  per  slot  are  independent  and 
identically  distributed  (i.i.d)  random  variables. 

We  define  the  delay,  experienced  by  the  n-th  arrived  packet,  as  the  time 

difference  between  its  arrival  at  the  transmitter,  and  the  end  of  its  successful 


transmission.  We  are  interested  in  evaluating  the  steady  state  statistics  of 
the  delay  process  {^n^n>l'  they  exist.  Due  to  space  limitations  in 

this  paper,  we  give  explicit  results,  only  for  the  first  moment  of  the  delay 
process.  However,  higher  moments  of  the  delay,  as  well  as  other  quantities  of 
interest  can  be  computed,  using  the  same  method* 

3.1  Example  1  ;  Controlled  ALOHA 

The  earliest  and  most  well  known  RMAAs  belong  to  the  class  of  the  ALOHA 
techniques  [13,6,16].  Here,  we  analyze  a  version  of  the  slotted  ALOHA 
algorithm,  that  operates  with  each  user  transmitting  a  newly  arrived  packet, 
in  the  first  slot  after  its  arrival.  Should  this  cause  a  collision,  each 
involved  user  independently  retransmits  its  packet  in  the  next  slot, 
with  probability  f. 

A  packet  whose  transmission  is  vmsuccessful  is  said  to  be  blocked.  Let  be 
the  number  of  blocked  packets  at  the  beginning  of  slot  1  (time  segment  [i,!-*-!)). 
This  number  will  be  referred  to  as  the  backlog  size.  Also,  let  denote  the 
number  of  blocked  packets  retransmitted  in  slot  i,  and  denote  the  number  of 
new  packets  transmitted  in  slot  1.  Given  >  m,  then  clearly, 

A 

P(Rj^-r)  -  BP(f)  -  (“)  (l-f)®">^  ,  i  -  0,1,2,...  (1) 

A  e"^  A" 

P(Ni-n)  -  Pn  -  -  .  i  -  0,1,2,...  (2) 

nl 

The  delay  process  induced  by  the  above  algorithm  "probabilistically 
restarts  itself"  at  the  beginning  of  each  slot  Tj^,  at  which  »  0,  i  «  1,2,...; 
this  is  so  because  the  number  of  arrivals  per  slot  is  an  i.i.d.  sequence  of 
random  variables.  Precisely,  let  T^  =  1,  and  define  T^^^^  as  the  first  slot 
after  T^  at  which  ^  =  0.  The  interval  (T^^,  i  =  1,2,...,  will  be 

referred  to  as  the  ith  session. 

Let  i  »  1,2,...,  denote  the  number  of  packets  successfully  transmitted 
in  the  interval  (0,  T^^,]  (  Mote  that  Rj^  also  represents  the  number  of  packets 
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arrived  during  the  interval  ^i  “  ^i+1  ”  ^  “  1,2,...,  is 

the  number  of  packets  successfully  transmitted  in  the  interval  (Tj^, 
sequence  is  a  renewal  process,  since  is  a  sequence 

of  nonnegative  i.i.d.  random  variables.  Furthermore,  the  delay  process 
is  regenerative  with  respect  to  the  renewal  process  (R *  with  regenera¬ 
tion  cycle, 

From  theorem  1,  with  we  have  that  if  C  =  e{C^}  <  •»,  and  if 

Cl 

S  B  e{  W^}  <  then,  there  exists  a -real  number  D,  such  that, 

iai 

In.  1  ^  .  S 

D  B  iim  -  tVj^  =  iim  -  E(  ZVj^)  »  e{  «}  =  -  a.e. 

n*«  **  i^l  n-»<"  i=1  ^ 

The  quantity  D  will  be  referred  to  as  the  mean  packet  de¬ 
lay.  Next,  we  develop  two  systems  of  equations,  whose  solution  may  be  used  to 
compute  the  mean  cycle  length  C,  and  the  mean  cumulative  delay  S*  The 
properties  and  the  computation  of  the  solution  will  be  postponed  until  section  4* 

l.a  Mean  Cycle  Length 

A  , 

If  the  mean  session  length  H  »  E  ^  ,  1  >  1,  Is  finite,  then 

by  Wald's  theorem  ,  we  have  that, 

C  B  IH  (3) 

To  determine  H,  we  proceed  as  follows.  Let  h^^  denote  the  random  number  of 
slots  needed  to  return  to  zero  backlog  size,  starting  from  a  slot  j  where  the 
backlog  size  is  equal  to  i,  i  >  0.  The  operation  of  the  alf^orithm  yields  the 
following  relation  for  the  h^'s. 


h 


0 


{ 


1 

J 


if  Nj  -  0,  1 
if  Nj  >  1 


(4. a) 
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+  h. 
1 


if  R.  +  N.  =  0 
J  J 


h.  =<  1  +  h_„  ,  I(i+N.-1>0)  if  R.  +  N.  =  1 

i+N^-1  j  j  3 


i>0 


1  +  h 


i+N. 
J 


if  R.  +  N.  >  1 
J  J 


(4.b) 


where  !(*)  is  the  indicator  function  of  the  event  in  the  parenthesis. 

If  we  let  =  E  {h^},  i  >  0,  then  after  taking  expectations  in  (4)  we  obtain, 

H  =  b,  +  E  c..  H.  >  i>0  (5) 

^  ^  k=0 


where  b^  =  1,  i^O; 

‘'oo  “  ‘^01  “  ‘^01  “  Pi‘  ‘^ik  =  Pk-i* 

‘'ii+1  “  =  PQ(1-Bjj(f))+  PjBj(f),  i>l; 

Cio  ■  0,  -  pQB^(f),  i>l;  =  0,  k<i-l,  i^l; 

where  p^^,  Bj(f),  ii^O,  p<j<i,  are  as  defined' in  (I),  (2),  respectively. 

Mote  chat  the  mean  session  length  H,  can  be  computed' from  system  (5),  since 
H  =  Hq. 


l.b  Mean  cumulative  Delay 

The  mean  cumulative  delay,  S  ,  can  be  computed  using  a  system  of  equations 
similar  to  system  ( 5 ) .  To  develop  such  a  system  we  proceed  as  follows •  Let 
denote  the  cumulative  delay  experienced  by  all  the  packets  that  were  successfully 

transmitted  during  the  h^  slots. Also,  let  •»  E{w^}j  i  ^0,  and  note  that 
S  =  Wq. 

The  operation  of  the  algorithm  yields  the  following  relations  for  the  w^'s. 


< 

Ni 

if 

N. 

=  0,  or  1 

w  «  < 

0  ( 

1 

if 

^3 

>  1 

(6. a) 


( 1) 


jUMLV 


Here,  for  convenience,  we  count  the  delay  of  a  packet,  starting  from 
the  beginning  of  the  first  slot  after  its  arrival. 


'•  'v'* 


i  +  N.  +  w. 
J  i 


if  R.  +  N.  =  0 
J  J 


w  =  /  i  +  N.  +  w, _  ,  I(i+N.-1>0)  if  R  +  N  =  1 


(6.b) 


i  +  N  +  w 

1  w  T 

.1 


if  R.  +  N.  >  1 
J  J 


After  talcing  expectations  in  (6)  we  obtain 


-  bj!  t  .  i>0 


where  b(|  =  X,bj[  =  i  +  X,  i>  1,  and  Cj^j^  are  as  defined  in  (5). 


3.2  Exaunple  2  :  The  "0.487"  Algorithm 


This  algorithm  is  the  most  efficient  RMAA  kno%m  to  date,  for  the  Poisson 


infinite-user  population  model  and  ternary  feedback;  (it  attains  a.maximum 


throughput  of  0.487  packets  per  slot) .  It  is  assumed  that  at  the  end  of 


each  slot  i  (time  segment  ,  the  us6rs --receive  feedback  “  0» 


or  c,  if  in  slot  i  there  were  respectively  zero,  one  or  more  than  one  packets 


transmitted.  For  the  description  of  the  algorithm,  motivation,  and  background 


discussions,  the  reader  is  referred  to  £81,  and  [15] 


Suppose  that  at  the  beginning  of  slot  v(time  segment  (v,v4-l]  ) ,  all  packets 


that  arrived  before  time  ty  <  v,  have  been  successfully  transmitted,  and  there 


is  no  information  concerning  the  packets  that  may  have  arrived  in  the  interval 


[t.^,v),  (i.e.,  the  distribution  of  the  interarrival  times  of  the  packets  in 


[tY,v)  is  the  same  as  the  one  ass\uaed  originally).  The  beginning  of  such  a 


slot  V  is  called  a  "collision  resolution  instant"(CRI) .  The  time  difference 


dy  =  V  -  ty  will  be  referred  to  as  the  "lag  at  v".  In  slot  v,  the  users 


that  generated  packets  in  the  interval  (t^,  t^  +  U^),  where  =  min(d^.  A), 


are  allowed  to  transmit;  A  is  a  parameter  to  be  properly  chosen  for  throughput 


maximization.  In  this  case,  we  say  that  the  interval  [tv«  +  Uy)  is  "trans¬ 


mitted".  After  a  random  number  of  slots  -f,  and  following  the  rules  of  the  algorithm. 


another  CRI  ,v',  is  reached,  with  a  corresponding  t^t>t  .  For  the  analysis  of  the 


algorithm,  we  need  the  following  definitions. 
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N 

ui 

V 

e{x|u} 


t 

V 

number  of  packets  in 

sum  of  delays  of  the  N  packets,  after  the  CRl  v 

sum  of  delays  of  the  N  packets,  until  the  instant  t^  +  U^. 

conditional  expectation  of  the  random  variable  X,  given  that 


Uv  =  u 


Let  be  the  sequence  of  successive  collision  resolution 

instants, and  let  dj^  be  the  lag  at  v^.  It  is  knovm,  [9  1,  that  the  sequence 

^^i^i>1  ^  Markov  chain,  with  state  space,  F  a  denumerable  dense  subset 

of  the  interval  (1,  “) .  Let  =  1,  d^  =  1,  and  define  »  as  the  first  slot 

after  T{,at  which  dm  -  1.  E^om  the  description  of  the  algorithm  it  can  be 

i+1 

seen,  after  a  little  thought,  that  the  induced  delay  process  probabilistically 
restarts  itself  at  the  beginning  of  each  slot  T^,  i  ■■  1,2, •••  .  Therefore, 
using  the  notation  and  definitions  of  example  1,  the  mean  packet  delay  D  is 
equal  to  S/C  provided  that  both  S  and  C  are  finite, 

2 •a  Mean  Cycle  Length 

As  in  example  1,  if  the  mean  session  length  H  =  e{T£^^  -  T^}  is 
finite,  then  C  =  To  evaluate  H  we  proceed  as  follows. 

Let  h^  denote  the  random  number  of  slots  needed  to  return  to  lag  equal 
to  one,  starting  from  a  collision  resolution  instant  v^  with  d^  x  d.  Note  that, 
by  definition,  h^  is  the  session  length.  The  operation  of  the  algorithm  yields 
the  following  relations  for  the  h^'s,  d  e  F. 

if  i  -  1 

(S.a) 

i+hd-6+)l  if  i  >  1 


d  >  A  ,  hjj  =  i+h^_5+g. 


(S.b) 


Taking  expectations  in  (8)  yields: 


1  <  d  <  a  , 


(9. a) 


5ii! 


i 

I 


=  E{A|d}  +  E  p(r,sld)Hd_r+s 


d  >  A  ,  Hjj  =  e{£|A}  +  Z  p(r,s|  A)Hjj_r+s 

s,r 


(9.b) 


where  p(r,s|x)  is  the  joint  conditional  probability  distribution  of  6,  and 
t,  at  the  point  values  r  and  s,  given  that  the  transmitted  interval  is  of 
length  x.  Note  that. 


p(r,l|x)  “  I 


(1+Xx)e"^*  if  r  =  X 


otherwise 


System  (9)  can  be  written  in  the  form 

Hd  -  bd  +  E  c^t  Ht  ,  deF  (10) 

teF 

where  b^  “  E{f|d},  1  <  d  <  A,  b^  =  E{f|A),  d  >  A,  and  where 

c^^,  d,  t  eF  are  nonnegative  coefficients  that  can  be  appropriately  identified 
from  (9).  nie  conditional  expectation  E{i|d},  1  <  d  <  A,  can  be  computed 
as  shown  in  Appendix  B. 

2.b  Mean  Oumulative  Delay 

Let  w^  denote  the  cumulative  delay  experienced  by  all  the  packets  that 
were  successfully  transmited  during  the  slots.  The.  operation  of  the  algorithm 
yields  the  following  relation  for  the  w^'s,  deF. 


1  <  d  <  A  , 


if  i  -  1 


“+<^d-6+i  if  *  >  1 


d  >  A  ,  “  •*>+<«•  (d-A)N+w^_g^^  if 


Taking  expectations,  we  obtain. 


'ri-Z'  V' 


11 


1  <  d  <  A,  w<j  =  E{u>ld}  +  E{<J»ld}  +  I  p(s,rld)W,j_j.+s  (11. a) 

r,s 

8^1 

d  >  A,  =  e{o)1a}  +e{«|»|A}  +(d-A)E{NlA}  +  E  p(s,r ( A)w^_r+3 

8,r 

(11. b) 

System  (11)  can  be  written  in  the  form 


■  b^  +  £  »  d  E  F  (12) 

•  te  F 

where  b^-stw  |d}+B{v|d),  1  <  d  <  A,  b4-E{«lA}+E{<»|A}  +(d-A)B{M| A)^ 
and  where  the  coefficients  d,  t  e  F  are  as  defined  in  (10).  The 

conditional  eiqMctations  B{tt|d}#  B(<|>|d}#^  1  <  d  <  A,  and  b{m|A}  can  be 
confuted  as  shorn  in  Appendix  B. 

4.  SYSTEM  SOLUTION  AND  MEAN  PAQCBT  DELAY  BOIIMDS 

In  this  section,  we  investigate  the  conditions  under  which  the  infinite 

dimensional  linear  systems..i54-f  ~(-7)^  4-iO  )-,~aad-(X2)::f  h«ve-4tnlqiie. -  . 

nonnegative  solutions,  and  we  develop  upper  and  lower  bounds  on  those  solutions. 
These  bounds  are  then  used  to  obtain  bounds  on  the  mean  packet  delay.  We  pro¬ 
ceed,  following  the  steps  outlined  in  section  2. 

4.1  Step  1 

EOr  convenience,  we  rewrite  an  infinite  linear  system  in  an  operator  form. 
Specifically,  let  E  be  the  space  of  sequences  X  «  {x(v)}:  A  -^R,  %diere  A 
is  a  countable  set.  Also,  let  be  the  subspace  of  E  for  which, 

E  |cu„x(v)(  <  •,  meA  ,  veA,  eR 
veA 

We  define  the  operator  L  {l^(x)}:  E^-^E,  as  follows. 

,  neA  ,  X  e  El*  .  bj;  e  R 

veA 


In  this  notation,  systems  (5),  (7),  (10),  and  (12)  can  be  written 
in  the  form, 

=  L(S^),  e  (13) 

We  are  interested  in  the  existence  and  uniqueness  of  nonnegative  points 
c  E^,  that  satisfy  (13  ) ;  such  points  will  be  referred  to  as  fixed 
points  of  L,  and  represent  solutions  to  the  corresponding  infinite  linear 
system  of  equations*  The  question  of  uniqueness  of  a  fixed  point  S^,  or 
equivalently  of  the  solution,  (s^(i)}#  to  the  system  that  operator  L 
represents,  depends  upon  what  conditions  are  imposed  on  the  solution*  Thus, 
after  the  existence  of  a  solution,  {s^(i)},  has  been  established,  one 
has  to  Indicate  a  class  of  sequences  in  which  the  solution  is  unique*  If  the 
algorithad.c  sequences  of  Interest  (Hj,},  or  (w^)  belong  to  the  indicated 
class,  then  they  must  coincide  with  the  solution  (s^d)}*  (This  will  be 
examined  in  Step  2)* 

Appendix  A  includes  a  number  of  results  that  can  be  used  to  establish 
existence  and  uniqueness  of  a' ^Ixed^^int  of  an  operator*  ~  Depend!!^  on  the 
operator,  some  are  more  straightforward  to  apply  than  others*  Among  the 
results  in  Appendix  A  that  can  be  used  to  establish  existence  of  a  solution. 
Lemma  A*  2  is  usually  the  most  useful*  According  to  Lemma  A* 2,  to  establish 
existence  of  a  nonnegative  fixed  point,  S^,of  a  nonnegative  operator,  L,  it 
suffices  to  find  a  point  e  E^,  such  that, 

0  <  L(X®>  <  X°  (14) 

A  point  X^,  satisfying  (L4  ) ,  also  serves  as  an  upper  bound  on  S^*  Further¬ 
more,  to  establish  a  lower  bound  on  S^,  it  suffices  to  find  a  point  e  E^, 
such  that, 

y°  <  L(Y°)  <  x°  (15) 

Thus,  under  (14)  and  (15),  we  have  that. 


L3 


Y®  <  <  X° 


(16) 


We  proceed  now  with  the  analysis  of  the  systems  developed  in  section  3. 

1.  Controlled  ALOHA 

System  (5)  —  Existence?  System  (5)  corresponds  to  an  operator  with 

Li  Li 

by  =  by,  Cyv  ■  Cyy,  U,  V  c  Nq,  whoro  is  the  set  of  nonnegative  integers, 

and  the  bj^’s  and  Cy^'s  are  as  defined  in  (5).  If  we  let  X°  ■  {x°(k)}  with  x°(0)‘=c^, 

x°(k)  =  ot^k  +  6^,  k  >1,  then  by  straightforward  manipulations  we  have  that,  for 

this  choice  of  X^,  (14)  is  satisfied  if  and  only  if  the  following  inequalities  are 

satisfied. 


1  <  €^(f)  -  every  k  >  1 

>  8up|— — -  ,  k  >  1  > 

Uk<«)-*  * 

6^  >  (1  -  o„  (Ci(f)-X))  /  (PqbJcO) 

<=u  >  1  +  %  <^-Pi)  - 


(17) 


(18) 


(19.  a) 
(19.b) 


It  can  be  readily  seen  from  (17)  that  if  the  retransmission  probability  f  is 
constant  in  every  slot,  then  there  is  no  A  >0  for  which  (17)  is  satisfied. 

If  the  retransmission  probability  f  ,  at  each  slot  i,  were  alloired  to  depend 
on  the  current  backlog  size,  in  accordance  to  a  stationary  control  policy 
f  «  f(M^),  then  it  is  of  interest  to  choose  f('  )  so  that  it  maximizes  the  set 
of  A's  for  which  inequality  (17)  is  satisfied.  This  is  equivalent  to 
maximizing  €)((f)  with  respect  to  f.  It  can  be  easily  verified  that,  for 
every  k  >1,  6jj(f)  is  maximized  for  f(k)  =  f*(k),  where 


2.  We  should  mention  that,  in  a  distributed  environment,  the  backlog  size 
dependent  retransmission  probability  f*(* )  is  nonimplemen table,  since 
users  are  not  aware  of  the  current  backlog  size.  However,  the  control 
policy  given  by  (20)  can  be  implemented  approximately  by  adaptive  control 
schemes  that  estimate  the  current  backlog  size  using  observable  feedback 
information  from  the  past  activity  on  the  channel  [6,16]. 


From  chls  point  on,  we  assune  that  C  is  chosen  as  in  (20)  •  Under  this  assumption 
inequality  (17)  is  satisfied,  provided  that. 


X  <  inf  {€jj(f*),  k  >  1}  -  e 


-1 


To  satisfy  inequalities  (18),^  and  (19),  we  choose. 


-  sup 


— — - .  k  >  1 


e-l-X 


(21. a) 


6^  -  e^  -a^  (1-X  eS  ,  -  1  +  X  (1-e^  )  +  6^^  d-e^  -  Xe^)  (21. b) 

Sifflildily,  it  is  straightforward  to  show  that  if  X  <  e’^,  then  the  point  Y°  with 


y^(0)  *•  c^,  y*’(k)  •  a^k  +  k  i  1,  and 

-  (^  e^-X  -  e^-  a^(l-Xe^),  -  1  +  o^X  (1-i^)  +  B^  (1-e^-Xe^) 

satisfies  (15).  Thus,  from  (16)  and  for  X  <  e  ^  we  have  that  system  (5)  has  a 

L  Lj^ 

solution,  S  >  (s  (k)},  such  that 


0  <  _<  s  ^(0)  1  c^;  0  <  a^k  +  B^  <  s  ^  (k)  <  a^k  +  B^  ,  k  k  1  (2! 

where  a^,  B^,  are  as  given  by  (21),  and  B^,  are  as  given  by  (22). 


^2 

System  (7)  —  Existence:  System  (7)  corresponds  to  an  operator  L2  with  b||  = 
L2 

Cpv  ^  M  ,  V  e  Kg,  where  the  b'^'s  and  Cp^'s  are  as  defined  in  (7). 


Due  to  the  fact  that  b]^  is  a  linear  function  of  k,  and  since 


L5 


it  can  be  easily  seen  that  there  is  no  linear  sequence =  {x°(k)}  satisfying 
(14).  However,  given  X  <  e  it  is  straightforward  to  show  that  we  can  choose 
coefficients  Y^»  such  that  the  point X°  with 

x°(0)  =  d^,  x°(k)  =  +  6^k  +  k  i  1,  and  the  point  y°(0)  =  d^, 

y°(k)  =  Y^k^  satisfy  (14),  and  (15),  respectively.  The 

following  is  such  a  choice: 

Y^  -  0.5  (i^-X)"^,  6^  -  2y^  (X+y^(X  +  X^  +  i^))  (24. a) 

+  Y^  (1  +  X  (1  +  X)  e^)  +  6^  (Xe^-  1)  (24. b) 

du  -  d  (Y^.  6^.  C^)  ^  X  +  X  (1+X  -  e"^)  Y„  +  X  (1-i^)  6^  +  (l-i^-Xi^)C^ 

(24. c) 

Y^  -  0.5  (^  S^-X)"^.  6^  -  2Y£  (X  +  Y^a  +  X^  +  i^(l-2X)))  (25. a) 

“  C  ^  (25. b) 

where  C(*.*)»  d(*,*,«)  are  as  defined  in  (24. b),  (24. c),  respectively. 

“1  -  ^2  ^2 
Thus,  if  X  <  e  ,  then  system  (7)  has  a  solution  S  ■  (s  (k)}  ,  such  that, 

0  <  d^  <  8  (0)  <  d^  (26. a) 

0  <  <  sM  <  Y^  +  «„k  +  ,  k  >  0  C26.b) 

where  Y^*  5^,  d^  are  as  given  by  (24),  and  yj^,  6^,  C^,  d^  are  as  given  by  (25). 

Syscems  (5)  and  (7)  —  Unicrueness 

^1 

We  will  show  that,  both  the  solution  (s  (i)}  of  system  (S)  and  the 
^2 

solution  (s  (i)}  of  system  (7)  are  unique  in  the  class 

|x(i)| 

X  :  sup  — —  <  »• 
iCNo  i^+c 

where  c  is  a  positive  constant. 

We  start  with  system  (7).  Since  L2  is  majorant  of  itself,  from  theorem  A.l, 


such  chat  0  <  S*  <  S 


we  have  that  L2  has  a  principal  fixed  point  S*  , 


According  to  theorem  A<2#  the  fixed  point  S  is  unicpie  in  the  class  / 


.^2  ( 

:*  -  j  X  : 


sup 

ICN. 


|x(i)| 


o  s^2(i) 


provided  that  e  .  since/ 


Y°  e  E2/  S  will  be  unique 
L2 

in  E2'  if  show  that  E*  *  According  to  lemma  A.  1,  it  suffices  to  show 

that/ 


L2 

8*  (1) 

sup  -- — —  <  •• 

ieM_  i2+c 


(27) 


and 


__SitJi)  - 

inf  - >  0 

iCN  i2+c 


(28) 


L2  1-2 


Since  0  <  a*  (i)  <8  (1)/  (27)  follows  from  (26 )  •  YO  show  that  (28) 

L2 

holds /  we  use  the  power  sequence/  (Sjj  ^2  *'i*^)*  initial  point  0.  By 


definition  (see  Appendix  A) /  Sn  is  the  point  that  results  after  operates  n 


1-2 


times  on  the  zero  point/  (i>e«/  S^  ■>  L2(0))/  and  S^^  S*  /  as  n  .  Due 
L2  1*2 

to  the  fact  that  >  0/  c^jj  >  0,  i,  k  e  Nq,  we  have  that  0  <  s,^  (i)  < 

s^+^d)  <  s*  (i)/  for  every  n  >  1/  i  >0.  Also,  it  can  be  readily  shown  by 


induction  that,  for  every  i  >1,  n  >1, 


s_  (i)  >  ni  -  n(n-1)/2 


(29) 


Fron  ( 29)  we  obtain 


L2  L2 

s*  (i)  Si  (i)  i2^i  1 

Aim  inf - - —  >  Aim  inf -  >  - 2 -  “  ~  (30) 

i-»<»  i2+c  i+*  i2+c  i"**®  2(1  +c)  2 

1*2 

(26)  follows  from  (30),  and  the  fact  that  s*  (i)  >  0,  i  >0. 

The  uniqueness  in  E2  of  the  solution  {s  (i)}  of  system  (5)  follows 
from  theorem  A.4,  part  (ii),  after  one  identifies  with  O2  and  L2  with  0^, 
in  the  theorem. 

2.  The  *’0.487*'  Algorithm 

System  (10)  —  Existence  and  Initial  Bounds 

Li  L, 

System  (10)  corresponds  to  an  operator  L^  with  by  ■  by  ,  Cy^  ■ 
ti,  V  e  F,  where  the  by-ls-and  Cy^'s  are  as  defined  in  (10).  To 

establish  t)ie_  existence  of  a  nonnyatlve  soluti«m-to  sy^^^  ***ll**"_ - 

the  same  procedure  as  in  system  ( 5 ) . 

Let  X°  -  {x°(d)}  with  x(d)  -  o^^d  +  8„,  d  eF  ,  and  let  X'  -  L,(X  ®). 

After  straightforward  manipulations,  we  obtain, 

X' (d)-x°(d)+E{A|d>+o^(E{Ald}-E{6|d}-(l+Ad)e"^‘^e^(l+Xd)e“^‘*  ,  1  1  d  <  A 

(31. a) 

x'(d)  -  xO(d)  +  E{A|d>  -  a^(E{6lA>  -  e{A1a>),  d  >  d  (31. b) 

According  to  Lemma  A. 2,  to  establish  the  existence  of  a  nonnegative  fixed  point 
of  L^,  it  suffices  to  show  that  there  exist  a^,  such  that. 


0  <  x'(d)  <  xO(d)  ,  for  every  d  e  F 


(32) 
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If  the  condition 


e{6|A}  >  E(i|A} 


holds,  then  it  can  be  readily  seen  from  (31)  that  (32)  is  satisfied,  if 


we  choose  a,,,  8^^  as  follows: 


eU|a) 


e{6|a}  -  E{f|A} 


(34. a) 


6^  »  max(-  o^,  sup  (p(d))} 
1<d<A 


(34. b) 


where 


ZU  |d}  +  a^(B{A(d}  -  E{6|d}  -  ( 1-»-Xd)e:9(-Xd) ) 


( 1+Xd)exp(-Xd) 


The  conditional  expectations  appearing  in  the  above  expressions  cem  be  computed 


as  shovm  in  Appendix  B« 


Similarly,  it  can  be  shown  that,  under  (33),  the  point  Y®  -  {y°(d)> 


with  yO(d)  =  o^^d  +  8^,  d  e  F  satisfies  the  inequality  Y®  <  I,(Y®)  <  X9, 


if  aj^  and  are  chosen  as  follows: 


»£  -  «u  ' 


inf  (  p(d)) 
1<d<A 


where  a^,  p(d)  are  as  given  by  ( 34) • 


Thus,  if  (33)  holds,  then  from  lemma  A. 2  we  have  that  system  (10)  has  a 
nonnegative  solution  S  ,  such  that. 


Ojjd  +  <s  (d)  <  Oyd  ♦  0^  ,  d  e  F 


where  a^,  and  Ojj,  are  as  given  by  (34)  and  (35),  respectively. 


System  (12)  —  Existence  and  Initial  Bounds 

Let  L2  be  the  operator  that  corresponds  to  system  (12).  Also,  let 


X®  -  {xO(d)}  with  x®(d)  -  Y^d^  +  6^d  +  d  eF  ,  and  -  {y®(d)}  with 


>(d)  =  Yffd2  +  6-d  +  ,  d  e  F. 


Ebllowing  the  same  procedure  as  for  system  (10),  we  can  show  that  if  ( 33 ) 


holds,  then  system  (12)  has  a  nonnegative  solution  S  >  (s  (d)},  d  e  F, 
such  that. 


Y^d^  -t-  6gd  +  <  a  (d)  <  Y„d^  ^u<* 


..  (373 . 


where. 


Bl^b}  __ 

2(e{6Ja}  -  eU|A}) 


e{o)|A}  +  e{4'|A}  -  A  e{h|a}  +  Y^  e{(6-)^|A} 
E{4|A}  -  E{t|A} 


=  sup  (♦(d))  ,  Ct  “  inf  (♦(d)) 
1<d<A  1<d<A 


♦  (d) 


2{u|d)+E{^|d}+Y^(E{(6-i)2|d}-2dE{6-i|d}-(l+Xd)e~^^)-6y(E{6-£,|d}+(l+Xd)( 

(1+Xd)c"^<i 


The  conditional  expectations  in  the  above  expressions  can  be  computed  as 


shown  in  Appendix  B. 


Remark  It  Is  known  .[7]  .that  Inequality  03  )  is  satisified  if  A  <  Aj^(A); 
where  is  maximized  for  A  «  2.6,  and  A^^^  (2.6)  =  0.4871. 


Systems  (10)  and  (12)  —  Uniqueness 

He  will  show  that  both  systems  (10)  and  (12)  have  unique  solutions  in  the 

class 


^2 


sup 

Atf 


(38) 


As  in  the  case  of  systems  (7)  and  (9)  in  example  1,  if  we  show  uniqueness  for 
system  (12),  then  the  uniqueness  for  system  (10)  follows  from  theorem  A.4, 
part  (ii)* 


According  to  fiheorem  A*  2, 


1*2 

the  fixed  point  S  is  unique  in  the  class 


1*2  (  l*<d)|  ) 

E*  -  <  X  :  sup  — r-—  < 

(  deF  s,i2(d)  ) 

1j2  I»2 

provided  that  V®  e  F*  .  Since,  by  construction,  Y®  ^  ^2'  ®  will  be 
1-2 

unique  in  E2,  if  E«  ^  ^2'  show  that  the  latter  holds,  we  proceed 

as  follows. 

Let  be  the  power  sequence  of  L2  witli  initial  point  0.  Clearly, 


sj^(d)  “  b^  >  e  >  0,  for  every  deF 


(39) 


Also,  it  can  be  easily  shown  by  induction  that. 


1*2  n(n“1 ) 

s„(d)>n((d-A)E(N|  A}+e{u|  A}+e{<|»|  A})-  ----- (e{  6 1  A}-e{  i  |  A)  )e{n|  A}  (^0) 


Etor  d  >  2A,  letting' 


for  every  d  e  F,  n  >  1 ,  such  that  d  >  nA  . 
d 


n=|^-j  -  t  in  (40)/  and  using  the  fact  that  [^  “  ^  yields. 


^  d 

5^(d)  >ad2  +  0d+Y/d>2A  ,n>[-J- 


(41) 


where  a  >  0.  (The  expressions  for  the  coefficients  u,  d,  Y  are  not  of 


interest  and,  therefore,  are  omitted) . 


If  S*  is  the  principal  solution  of  L2,  then  from  lemma  A«2  we  have 


s*  (d)  >  8j^(d)  >  0  ,  for  every  d  e  F  ,  n  >1 


From  (  39)  and  (42)  we  have  that. 


2  2 

s*  (d)  >  max  (e,  ad  +  3d  +  y)  ,  V  dcF 


From 


(43;  we  conclude  that. 


^2 

s*  (d) 

inf  - >  0 

de  F  d2 


From  (37),  and  the  fact  that  S*  <  S  ,  we  have. 


L2 

s*  (d) 

sup  - <  ® 


deF 


Finally,  from  (44),  (45),  and  lemma  A. 1  we  have  that 


E.  =  E. 


3.  laj  denotes  the  maximum  integer  not  exceeding  a. 


(42) 


(43) 


(44) 


(45) 


4 . 2  Step  2 


In  step  1,  we  have  established  conditions  for  the  existence  of  nonnegative 
solutions  to  the  systems  of  interest,  and  we  have  identified  classes  of 


sequences  in  which  these  solutions  are  unique.  Here,  we  show  that  the  algorithmic 
sequences  where  Hj_  =  E  {hjj^}  and  =  E  {wj_} ,  belong  to  the 

corresponding  identified  class,  and  therefore,  coincide  with  the  unique  solution 
in  the  class.  The  proof  is  based  on  theorem  A. 6,  and  is  the  s^une  for  the  two 
algorithms. 

Etor  the  case  of  the  sequence  let,  in  theorem  A. 6,  L  ®  L'j, 

Xf  =  h^,  and  Xj  =  min(h^,  n),  n  =  1,2,3,...  .  By  definition,  the  X^*s  and 
xj's  satisfy  condition  (a)  in  the  theorem.  Condition  (b)  follows  from  the 
fact  that  x"  <  n  a.e.  .  Finally,  condition  (c)  follows  from  the  operation 


of  the  algorithm.  Thus,  {u^}  ~  S 

L2 

Similarly,  to  show  that  {Hj^}  <=8  , 

Xj^  =  Wj^,  and  x"  “  mln(w^,  n) ,  n  =  1,2,3, 


we  apply  theorem  A* 6, 

«  «  *  » 


with  L  L2, 


4.3  Step  3 

o  o 

In  step  1,  we  have  already  found  upper  and  lower  bounds,  X  emd  Y  , 
respectively,  on  the  solutions  to  the  systems  of  interest.  These  bounds 
can  be  improved  either  by  computing  the  power  sequences  of  the  corresponding 
operators  with  initial  points  the  bounds  X°  and  Y°,  (lemma  A. 2),  or  by  solving 
finite  systems  of  linear  equations  that  are  truncations  of  the  original  in¬ 
finite  systems,  (theorem  A. 5).  Both  methods  can  provide  arbitrarily  tight 
upper  and  lower  bounds.  We  use  the  first  method  in  the  "0.487"  algorithm, 
and  the  second  method  in  the  controlled  A1/}HA. 


1 .  Controlled  ALOHA 

For  system  (5),  we  apply  theorem  A. 5  with  L  =  and. 


i 


U  (1)  =  Oyl  +  Uu 


,  i  e  N, 


I  (i)  =  ,  i  e  N^. 


A.  =  {0,1,2...,j}  /  j  e  Nq 


;where  ot^^,  and  are  as  given  by  (21)  and  (22),  respectively.  Note 

that,  for  j  <  <*>  ,  Aj  is  a  finite  set  and,  therefore,  all  conditions  in  the  theorem 


are  satisfied.  Thus,  for  X  <  e”^. 


s  (i)  <  Hi  =  s  (i)  <  s  (i)  ,  0  <  i  <  j 


where  {s  and  {s  (i)}o<i<j  unique  solutions  of 

the  (j+1 ) -dimensional  systems  (46)  and  (47),  resjpectively. 


Gj  j 

Hi  “  bi  +  E  Cij^  ,  i  0  <  i  <  j 


hJ  =  bi  +  E  Ci,^  H*  ,  i  0  <  i  <  j 
k=0 


0^ 


where  b^  ,  b^  are  as  defined  in  the  theorem  with  Pi=ai=bi,0<i<j. 


therein.  We  solved  systems  (46)  and  (47)  for  j  =  50.  The  obtained  upper  bound 


u  S, 

Hq  and  lower  bound  on  the  mean  session  length  H^,  can  be  found  in  table  1, 
for  different  values  of  A,  (A<e^).  Bbr  system  (7)  we  followed  the  procedure 


described  above  with. 


h  =  hr 


u  (i)  =  Yu  i  +  6u  ^  +  '^u  '  i  ^  Ho 


^  (i)  =  Yi  i  +  <5^  i  +  ,  i  €  Nq 


Pi  »  Oi  =  bi  ,  i  e  Aj  =  {0,1,2,...,j}  ,  j  e  N<. 


kV  d  r  «  * 


t 


where  Yu»  Cu  are  as  given  by  (24),  and  Yj^/  Cjj  are  as  given  by  (25)  •  The 

obtained  bounds  W^,  on  the  mean  cumulative  dealy  W  are  included  in  table  1; 
o  o  o 

they  were  computed  using  j  SO.  ETom  the  regeneration  theorem  and  (3)  we  have^^ 


D  - - +  0.5 


The  upper  bound  D'*^  =  +  0.5,  and  the  lower  bound  D*'  •*  W^/{Xh“)+0.5 

on  D  are  included  in  table  1.  .Note.that.«  according  iu)  theorem  h.5,  arbitraurily 
tight  bounds  can  be  obtained  by  increasing  j.  ETom  a  theoretical  view  point  the 
bounds  become  exact  as  j  -*-  **. 


2.  The  *0.487**  Algorithm 

^1  _ 

Et-om  section  4.2  %ra  have  that,  £or  X  <  0.487,  s  (d),  d  e  r,  and 

L2  h-y  L2 

>  s  (d),  d  e  r,  where  S  and  S  are  the  fixed  points  identified  in 
section  4.1.  According  to  lemma  A. 2  we  have  that, 

l"(Y°)  <  S  <  L!f(X^)  ,  n=1,2,...,  d  e  F  (49) 

^2 

L§(Y2°)  <  S  <  I'2^X2*  '  n=1,2,...,  d  e  F  (50) 

»h.r.  JCO  -  (.„a  .  V  ,  e  F-  *1°  -  (V  *  d  e  F 

*1  -  t’fu  a  e  F-  *1  ■  d  c  F 

and  where  o^,  Y^,  Qg,,  Y,,,  6^^,  are  as  given  by  (34),  (35),  and 

(37).  Ebr  n  1,  and  d  -  1,  (49)  yields  the  following  bounds  on  the  mean  session 


length 


4.  The  additional  0.5  units  of  time  represent  the  mean  delay  of  a  packet,  until 
the  beginning  of  the  first  slot  following  its  arrival.  (See  footnote  1). 


where 


h"-b{*|i}  +  o^j(l-(l+A)e"^  +E{Jl|l}-E{5|l})+|l^^(l-(l+A)e  "^) 
»  h’J  -  (6^  -  ejt)(1-(l+X)e“^) 


The  above  bounds  can  be  found  in  tiable  2,  for  different  values  of  X,  (  X  <  0.487). 
Hsr  n  »  1,  and  d  »  1,  (50)  yields  the  following  bounds  on  the  mean  cumulative  delay 
over  a  session 


where 

w!J«-e{o»|  1)+b{<(>|  l)+Y„(l-(  1+X)e"VE{(f-iQ?|  1}-2e(6-X|  i}) 
+  fiy^(l-(l+X)e“*  -E{«-X|l})+c„(l"(l+X)e“^) 

wj  -  Wl  -(C„-Cji)(l-(1+X)e"^) 


nie  bounds  w!|  and  are  include<i'l:n  table  2.  From  the  regeneration  theorem 
we  have  D  -  Wj^/(Xh^).  The  upper  bound  o'*  »  w“/(XH^^),  and  the  lower 
bound  D*  -  /(Xh’J)  on  the  mean  packet  delay  D  are  included  in  table  2.  The 

upper  bound  Is  plotted  in  figure  1,  together  with  the  same  bound  for  the  contro¬ 
lled  ALOHA.  We  note  that  tighter  bounds  can  be  obtained  either  by  evaluating  the 
bounds  given  by  (49)  and  (SO)  for  higher  values  of  n,  or  by  the  method  of  truncated 
systems  used  in  the  previous  ex2unple.  In  both  methods,  however,  we  must  first 
compute  the  conditional  probabilities  p(4,f|x)  defined  in  (9),  which  is  a 
computationally  con^lex  task.  Note  that  for  the  found  bounds,  (i.e.,  for  n  »  1  in 
( 49)  and  ( 50) ) ,  such  a  computation  is  not  required. 


5.  CDHCLUSIONS  AND  PRIOR  WORK 

In  this  paper  we  have  Introduced  a  method  for  the  delay  analysis  of  RMAAs,  in 
which  the  induced  packet  delay  process  is  regenerative,  and  we  have  demonstrated 
its  wide  applicability  by  applying  it  to  two  specific  examples.  Ihe  method  is 
based  on  a  well  known  result  from  the  theory  of  regenerative  processes,  which  relates 
the  asymptotic  statistics  of  such  processes  to  quantities  that  refer  only  to  one 
cycle  of  the  process.  The  per  cycle  quantities,  (e.g.,  mean  cycle  length,  expectation 
of  the  sura  of  the  values  of  the  process  over  a  cycle) ,  are  evaluated  from  the 
solution  of  infinite  dimensional  systems  of  linear  equations. 

In  applying  the  method  to  the  two  ex2unple-algorithms ,  we  have  put  emphasis 
on  the  methodology  and  rigorous  derivations  rather  than  finding  short  cuts  in  the 
analysis  of  a  particular  algorithm.  In  doing  so,  the  essential  simplicity  of  the 
method  may  have  been  obscured.  However,  to  appreciate  the  simplicity  of  the  method, 
we  note  that  only  by  using  Lemma  A. 2,  one  can  obtain  with  minimal  effort: 

1)  A  lower  bound  on  the  maximum  input  rate  that  an  algorithm  maintains  with 

finite  delay,  (i.e«,  a  lower.. bound -oo-jthe-  maximum  stable- throughput-induced 
by  the  algorithm).  Note  iaiaf^r  th^cW  of  this  j^per,  the 

found  bound  coincides  with  the  maximum  stable  throughput* 

2)  Optimal  algorithmic  p«amet«P?^iMwraf^l^T^  retransmissIonSprebS^-^^^-”^^ 
bility  policy  in  the  ALOHA  algorithm,  and  the  -window  size  A  in  the  "0.487" 
algorithm) . 

3)  Initial  bounds  on  the  mean  packet  delay,  that  can  be  used  (if  so  desired) 
to  form  finite  linear  systems,  whose  solution  can  yield  arbitrarily  tight 
bounds  on  the  mean  packet  delay* 

The  algorithms  that  served  as  examples  in  this  paper  have  been  analyzed  in  a 
number  of  studies.  ETom  the  literature  on  AI/>HA-type  algorithms,  we  mention  the 
work  in  (6] ,  where  the  stability  properties  of  the  version  of  the  controlled  AIAHA 
algorithm  considered  here  have  been  studied,  using  a  Markovian  model.  Bie  optimal 
retransmission  policy  was  derived  in  (6)  using  Fake's  lemma,  but  the  delay  analysis 
problem  was  not  addressed. 

The  delay  characteristics  of  the  "0.487"  algorithm  have  been  studied  in  [9  1 , 
using  a  different  approach.  In  contrast  to  the  method  in  [9  ] ,  the  method  proposed 
here  does  not  require  the  computation  of  steady-state  probabilities  of  the  under- 
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lying  Markov  chain  and,  therefore,  it  is  computationally  simpler.  Furthermore, 
since  our  approach  is  based  on  the  asymptotic  properties  of-  regenerative  processes, 
it  yields  stronger  convergence  results. 


X 

”o 

**0 

”o 

0 

0.05 

1.00426 

1.00426 

0.05782 

0.05782 

1.65163 

1.65163 

0.10 

1.01983 

1.01983 

0.14067 

0.14067 

1.87936 

1.87936 

0.15 

1.05361 

1.05361 

0.27541 

0.27541 

2.24265 

2.24265 

0.20 

1.12017 

1.12017 

0.53225 

0.53225 

2.87576 

2.87576 

0.25 

1.25676 

1.25676 

1.14710 

1.14710 

4.15097 

4.15097 

0.30 

1.59883 

1.59883 

3.39345 

3.39345 

7.57485 

7.57485 

0.35 

3.48032 

3.51077 

37.04013 

39.39037 

30.64403 

32.83714 

Table  1 

Delays  for  bhe  OOnrrolled  ALOHA 


X 

L 

“i 

t 

w 

W« 

1 

0 

D" 

.01 

1.00025 

1.0003 

.015258 

.015258 

1.5253 

1.5255 

.05 

1.00395 

1.00474 

.08234 

.082346 

1.6348 

1.6388 

.1 

1.025 

1.030 

.18503 

.1859 

1.796 

1.8130 

.15 

1.060 

1.061 

.3212 

.3245 

2.040 

.2 

1.1167 

1.11367 

.5162 

.5254 

2.270 

2.352 

.25 

1.2069 

1.240 

.8243 

.8468 

2.66 

2.80 

.3 

1.356 

1.408 

1.381 

1.434 

3.525 

.35 

1.627 

1.710 

2.6088 

2.7423 

4.358 

4.8151 

.40 

2.2279 

2.374 

6.64  38 

6.8603 

6.779 

7.670 

.45 

4.487 

4.8536 

35.012 

37.871 

16.030 

18.754 

.47 

9.110 

9.916 

163.698 

178. 178 

35.125 

41.613 

00 

• 

21.175 

23.122 

944.35 

1031.12 

85 . 086 

101.452 

Table  2 

Delays  for  the  "0.487"  Algorithm 


APPENDIX  A 


We  present,  in  a  generalized  format,  some  basic  results  regarding  the 
approximate  computation  of  solutions  of  infinite  dimensionality  linear  systems 
[17].  Let  A  be  a  denumerable  set  of  indices,  and  let  E  be  the  space  of  sequences 
X  =  {x(k)}:  A  R.  Given  a  set  {  bj[|  CR,i,kcA  },  let  E^  be  the  subspace 

of  E  defined  as  follows:  E^  =  {X:  £  |c?'.x(k)|<  “,VieA  }.  we  define  an  operator 

keA 

L  :  E^  E  as  follows:  y(i)  =  Lj^(X)  =  bj^  +  £  c^^x{k),i  eA  ,xeE^.  A  point 

keA 

gL  egl*,  such  that. 


=  L(S^) 


(A.1) 


Is  called  a  fixed  point  of  the  operator  L*  (A.1)  represents  an  infinite  system  of 

linear  equations  and  a  fixed  point  is  a  solution  to  this  system.  Given  an  operator 

e  A  ^  n 

L,  we  define  its  n-th  power  l”  as  follows:  l'(X^)  ■  L(Xq)/  L  (X^)  ■  Jj(h  (X^))# 

n-1,2,...,  provided  that  X^  e?*,  and  l”(X^)  e^,  for  every  n  >  1.  The  sequence 
A  _ 

{x,j}  “  {l  (Xq)},  n»1,2,...  is  called  the  power  sequence  of  L,  with  initial 
point  Xq.  a  fixed  point  of  L  that  is  a  pointwise  limit  of  the  power  sequence  of  L, 
with  initial  point  Xq  =  0,  is  called  a  principal  fixed  point  of  !»,  and  is  denoted  by 
S«.  An  operator  0:  E^  >E  is  called  a  ma jorant  of  L,  iff. 


<  b 


0 

'ik 

0 


i,k  e  A 
i  e  A 


In  this  case,  L  is  called  a  minorant  of  0«  The  notation  X  <  X' ,  X  <  X' ,  X,X*eE 
means  that  x(k)  <  x'(k),  {x(k)  <  x*(k)),  k  eA  .  A  point  X  eE  is  called  positive 
(nonneqative)  iff,  0  <  X  (0  <  X).  By  jx]  we  denote  the  sequence  defined  by 


|x|(k)  a  |x(k)|,  k  c  A.  Theorems  Al,  A2  below  are  essentially  theorems  I,  II  ^  2 


of  [^7].  They  relate  the  existence  and  uniqueness  of  a  fixed  point  of  L,  to  the 


existence  of  a  fixed  point  of  a  majorant  6  of  L. 


Theorem  A . 1  If  G  is  a  majorant  of  L,  and  0  has  a  nonnegative  fixed  point  S  , 


then  both  0  and  L  have  principal  fixed  points  s|^«  Moreover,  0  <  |s^|  < 


Theorem  A. 2  If  0  is  a  majorant  of  L,  and  0  has  a  nonnegative  fixed  point 


0  Xj 

then  the  principal  fixed  point  of  b  is  unique  in  the  class  E*  CE  t  defined  as 


follows. 


E®  »  {x  eE  ;  sup  <  •}  1 1 ) 

ieA  s*(i) 


EUrthemore,  S*  is  the  pointwise  limit  of  £uiy  power  sequence  of  L,  with  Initial  point 
any  point  in  E^. 


Theorem  A. 3  below  relates  the  existence  and  uniqueness  of  a  fixed  point  of  L, 


to  the  existence  of  a  fixed  point  of  a  majorant  @  of  L,  and  it  is  a  consequence  of  the 


theory  of  regular  systems  [17].  Its  difference  from  theorems  A1,  A2,  lies  in  the  fact 


0  0 

that,  under  the  stated  assumptions  in  it,  we  have,  S  ■=  S*. 


Theorem  A. 3  If  0  is  a  majorant  of  L,  and  0  has  a  positive  fixed  point 


such  that. 


inf  - - >  0, 

ieA  s  (i) 


then  S*  =  S®.  Therefore,  theorem  A. 2  holds  with  S®  replaced  by  S®. 


( 1 )  We  adopt  the  convention : 


0  a 

-Ml,  -  =  a  >  0 


.3 


The  following  theorem  relates  the  existence  and  uniqueness  of  a  fixed  point 
of  some  operator  O2  to  the  existence  and  uniqueness  of  such  a  point  for  another 
operator  0^ r  where  the  latter  is  not  necessarily  a  majorant  of  the  former. 

Theorem  A. 4  Let  0^,  O2,  be  two  operators  such  that, 

0  0 

12 

(a)  Ci,^  >  |ci)cl  V  i,k  eA  ,  bi  e  lO,“),Vi  eA 


(i)  If  Of  has  a  fixed  point  S  , -and  there  exists  a  sequence  qrA-*-  R,  such  that, 

(b)  g  +  S  >0 

(c)  I  g{k)  I  <  «,Vi  e  A 

keA 

O2  Of  Of 

(d)  (b^  I  <  (bj|_  +  g(i)  -  £  g(k))  M,  VI  eA  ,  for  some  M  >  0. 

kcA 

O2 

Chen,02  has  a  fixed  point,  S  .  ^ 

(ii)  If  (a),  (b),  (d)  hold,  for  g  -  0,  then  S  is  unique  in  the 

Of  Of 

class  E*  ,  where  E*  is  as  defined  in  Th.  A. 2. 

(iii)  If  in  addition  to  (a),  (b),  (d),  we  have  that. 


'1 


0i  0, 

bj^  +g(i)-  £  Cij^g(k) 

keA 


(e)g+S  >0,  and  inf  - - - >0, 

ieA 

3  {i)+g(i) 

O2 

then  the  fixed  point  S  of  O2  is  unique  in  the  class  Eg  d E  ,  defined  as 


0,  O2 


follows. 


{xeE  :  sup 
ieA 


lx(i)| 


(? 


<  ®} 


s  (i)+g(i) 


A 

S  is  the  pointwise  limit  of  any  power  sequence  of  O2,  with  initial  point  in  Eg 


A.4 


Proof 

0i  C»1  0i  0i 

Part  (i):Let  Y  =  (S  +g)  M.  Since  S  =  ),  we  have  that, 

(?1 

Y  /M  -  g  =  Oj  (Y  /M  -  g) 
or 

0^  0\  0^  0^ 
y  (i)  =  M(bj^  +g(i)  -  ^  c^j^gCk))  +  Z  y  (i)  (A. 2) 

keA  keA 

ETom  (A. 2)  and  (b),  we  see  that  the  operator  0  with  par2uneters, 

0\  0l 

b®  =  M  (b^  +g(i)-  £  g(k)),  i  eA 

keA 

<=  ik  “  Cik  ^  ' 

has  a  nonnegative  fixed  point  S®  >■  Y  •  Because  of  (a)  and  (d),  0  is  a  Mjorant 

of  02>  f^om  theorem  A.1,  we  conclude  that  O2  has  a  fixed  point* 

0  0, 

Part  (ii)  :  This  follows  from  theorem  A. 2,  by  observing  that  S^=  and 

0  0-, 

therefore, 

Part  ( iii )  :  under  condition  (e),  theorem  A* 3  is  applicable,  and  shows  the  uniqueness 
of  the  fixed  point  in  Eg  • 

The  following  lemma  is  useful  in  identifying  the  class  within  which  the  fixed 
point  of  an  operator  is  unique,  in  the  case  where  the  solution  of  the  majorant  is  not 
exactly  known* 

Lemma  A*1  If  {s(l)}  ,  {({(i)}  :  A  -*■  R,  and, 

(®)  {s(i)},  {({(i)}  are  nonnegative 


s(i) 

(b)  sup  - <  «" 

ieA  ^(i) 


8(i) 

(c)  inf-—- —  >  0, 
ieA  6(i) 


lx(i) 1 

then  sup - <  " 

ieA 


lx(i) 1 

iff  sup - <  “ 

ieA  s(i) 


X  (£  ,  i.e.  the  classes 


r  r  P  ,  r  f  P  , 

Cs  “  Ixec:  sup - <  *}  and  =  Ixet*  sup  — - - <  «}  coincide. 

ieA  s(i)  ^  ieA  ({(i) 

Proof  For  the  "if"  part  let 
|x{i)| 

sup  -  =A<  (x(i)|  <As(i),  i  eA 

ieA  s(i) 


(A. 3) 


Because  of  ( b ) ,  we  have , 


s{i)  <  B  ({(i),  i  eA  ,  B  <  «  (A. 4) 

•  I  .  1  |x(i) I 

From  (A. 3),  (A. 4),  we  conclude  that  |x{i)|  <  A  B  o(i),  i  eA  ,  or  sup - <  A  B  < 

ieA 

The  proof  of  the  "only  if"  part  is  similar. 

The  lemma  below  is  used  to  establish  the  existence  of  a  fixed  point  of  an 
operator  L,  as  well  as  upper  and  lower  bounds  on  S^.  Monotonicity  is  proved  by  indue 
tion,  while  the  existence  of  a  fixed  point  is  established  via  the  extended  monotone 
convergence  theorem. 

Lemma  A. 2  Let  L  be  an  operator  with  nonnegative  parameters*  i.e.  : 

L  L 

‘^ik  ^  ®  i,k  eA  ,  bj^  >  0,  ieA,  if  there  exist  points  Y°,  X°  eE^, 

such  that. 


(a) 

Y° 

<  X° 

(b) 

X° 

>  L(X°)  >  0 

(c) 

Y° 

<  L{Y°), 

then  the  power  sequence  of  L,  with  initial  points  X°  (Y^) , decreases  (increases) 
monotonically  and  pointwise,  to  a  fixed  point  S^(S^).  Furthermore,  Y°  <  <  X°. 

and  ^  0. 

It  is  generally  difficult  to  establish  tight  bounds  on  S^,  using  the  method 
exhibited  by  lemma  A. 2.  The  following  theorem  provides  an  alternative  method  for 
the  computation  of  such  bounds.  Its  proof  is  based  on  theorems  A.l,  A. 2,  and  A. 3, 
and  is  straightforward. 
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(d)  If  in  addition  to  the  previous  conditions,  also  S  >0,  and 
‘^i 

inf  -  >  0,  then  the  operators  ♦j,  Fj,  0j  have  respective 

ieA  s^(i) 

^  ♦j  Fj 

unique  fixed  points,  S  ,  S  ,  S  »  in  the  class 


E  =  <  X  tE  ;  sup 

ieA 


x(i)|  ) 

- <  ooi  , 

^3.  .  I 


and  S 


Fj  Oj 

<  S  <  S 


s  (i) 


Remark  If  A^  is  a  finite  set  with  >  0,  Vi  e  A j ,  the  conditions  In  (c)  and  (d)  are 
clearly  satisfied.  If  in  addition,  ~  A^^A,  then  it  can  be 


L  L 

shown  that,  S*  -*•  S  ,  and  S*  -*■  S  ,  pointwise. 


The  quantities  of  interest  in  the  various  random  access  algorithms  are  statistics 
of  random  variables,  where  mauiy  of  those  statistics  are  fixed  points  of  some  operator  L. 
Theorem  A. 6  is  used  to  justify  the  latter  statement  and  appeared  in  [14], 

Theorem  A. 6  Let  L  be  an  operator  with  nonnegative  pareuaeters  that  has  a  unique 


L  L  ^ ' 

nonnegative  fixed  point  S  in  the  class  Eg  =  {xeE  :  sup  — —  <  “>} . 

■  ieA  g(i) 


Let  txj^*^},  txj^},  i  cA  ,  n  e  N,  be  families  of  random  variables,  such  that, 

(a)  0  <  x?^x^,  a.e.  for  every  i  eA 

(b)  X?  <  M^g(i),  a.e.  for  every  i  eA  ,  <  “ 

(c)  f  <  L(f  ),  f  =  L(f),  where  f  (i)  =  E  (x. ,  f(i)  =  E  {x. } 

n  n  n  1  2. 

Then,  f  coincides  with  the  unique  fixed  point  in  E^. 

Proof 

We  observe  that  because  of  (b),  f  eE^,  and  because  of  (c)  and  lemma  A. 2,  f  <  S^. 

n  g  ’  n  — 

From  (a)  and  the  monotone  convergence  theorem,  we  have  that  f  increases  to  f 

n 

pointwise;  thus,  f  <  S^,  which  implies  that  feE^.  The  assertion  now  follows 

“  g 


from  the  fact  that  f  is  a  fixed  point  of  L. 


B.l 


APPENDIX  B 


In  section  4.2,  we  saw  that  the  computation  of  conditional  expectations, 
E{x|d},  is  required.  In  this  appendix,  we  show  that  those  conditional  expectations 
can  be  computed  with  high  accuracy.  Let  us  define. 


E{xld,  k}  :  The  conditional  expectation  of  the  random  variable  X,  given  that 
the  arrival  interval  contains  k  packets,  and  has  length  d. 


Then, 


E{X|d}  =  Z  E{xld,k}e“^‘^ 
k=0 


(B.l) 


Using  the  rules  of  the  algorithm,  the  quantities  E{X|d,k}  can  be  computed 
recursively  ,  as  follows. 


E{il/d,k} 

E{)1/1,0} 

E{)l/l,k} 


E{il/l,k};  V  deF 
E{i/l,l}  »  1 


(B.2) 


k-1 


(l+pf+E{)l/l,k-l}Pi+  Z  E{it/l,i>pJ)/(l-2pJ);  k  >  2 

i=2 

where  P^  =  (^)2  ^ 


E{6/d,k} 

E{6/1,0} 

E{6/l,k} 

E{Jl^/d,k} 


d  E{6/l,k};  V  deF 
E{5/1,1}  =  1 


k-1 


(B.3) 


(P^+P^:E{6/l,k-l}Pi+  z  E{6/l,i}P^)/(2-(l-p5));  k  >  2 
^  i=2 


E{il2/i,k};v  deF 
E{i^/1,0}=  =  1 


E{Jl  /l,k}= 


k-1 


(2E{Jl/l,k}+2E{ll/k,k-l}*p!^+pl!^-l+pl['E{i^/l,k-l}+  Z  E{jl^/l,i}pb/ 

^  ^  ^  i=2  ^ 

(l-2Po);k>2  (B.4) 


E{5^/d,k}=  d^E{6^/l,k}  ;  v  deF 
£{6^/1, 0}=  E{62/1,1}  =  1 

E{6^/l,k}=  (.25(P^Pq)+.5  PQE{6/l,k}+.5  P^  E{6/l,k-l}+.25  E{6^/l,k-l}P^+ 


k-1 


.25  Z  E{6^/17i}P^)/(l-.5pI^)  ;  k  >  2 
i=2  ^  ^ 


(B.5) 


E{6£/d,k}  =  d  E{61l/l,k}  ;  V  deF 
E{6)1/1,0}  =  E{6JI/1,1}  =  1 

E{6Jl/l.k}  =  (E{6/l,k}+.5pjfE{£/l,k-l>+.5P^E{Jl/l,k  }+.5  P^+.5E{6/l,k-l}P^+ 

+.5P^E{6£/l,k-l}+.5  E  E{5Jl/l,i}p5^)/(l-ph  ;  k  >  2  (B.6) 

i»2  ^  ° 

E{N/d,k}  =  E{N/l,k};  V  deF 
E{N/1,0}  =  0,  E{N/1,1}  =  1 

E{N/l,k}  =  P*^+P^{N/l,k-l}+*^E^E{N/l,i}p^;  k  >  2  (B.7) 

i=2  ^ 

EWd.k}  »  E{a)/l,k};  V  deF 

E{cj/1,0}  «  0,  E{w/l,l}  -  1  (B.8) 

k— 1 

E{w/l,k}  -  (P^+E{N/l,k}+P^E{N/l,k-l}+p!^{a)/l,k-l}+  Z  E{w/l,k}ph/(l-2pb  ;  k  >  2 

i_2  I  u 

E{^()/d,k}  =  d  E^/l.k};  V  deF 

E4/1,0}  «•  0,  E{i|>/l,l}  -  j  (B.9) 

^  k-1 

-  (E{N/l,k}(l-pJ)-pjE{N/l.k-l}+.5P^4E{4>/l,k-l>P^+  S  E{tf»/l,l}pb 

k  ^  i-2  ^ 

/(2  (I-Pj))  ;  k  >  2 

From  formulas  (B.2)-(B.9)y  we  see  that  a  finite  number,  M,  of  terms  from  the 

infinite  "series "(B. 10 cam  be  easily  computed.  Also,  for  large  k 'values,  and  based 

on  the  recursive  expressions,  simple  upper  and  lower  bounds  on  E{X/d,k}  can  be 

00  —Id 

developed.  Those  bounds  can  be  used  to  tightly  bound  the  sum  E  E{X/d,k}e 

k-M+1 

Remark  It  can  be  also  proved  that 
E{N/d}  -  XE{6/d} 

E{i<>/d}  -  Xd  E{6/d}  -  XE{6^/d)0.5 
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